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We construct a one-dimensional totally asymmetric simple exclusion process (TASEP) on a ring 
with two segments having unequal hopping rates, coupled to particle non-conserving Langmuir 
kinetics (LK) characterized by equal attachment and detachment rates. In the steady state, in the 
limit of competing LK and TASEP, the model is always found in states of phase coexistence. We 
uncover a nonequilibrium phase transition between a three-phase and a two-phase coexistence in 
the faster segment, controlled by the underlying inhomogeneity configurations and LK. The model 
is always found to be half-filled on average in the steady state, regardless of the hopping rates and 
the attachment/detachment rate. 


I. INTRODUCTION 

Totally asymmetric simple exclusion process (TASEP) 
and its variants with open boundaries in one dimension 
(ID) serve as simple models of restricted ID transport. 
These ID transports are observed in a variety of situa¬ 
tions, e.g., motion in nuclear pore complex of cells [lj, 
motion of molecular motors along microtubules Q, fluid 
flow in artificial crystalline zeolites Q and protein syn¬ 
thesis by messenger RNA (mRNA) ribosome complex in 
cells [H; see, Refs. Q for basic reviews on asymmetric ex¬ 
clusion processes. The coupled dynamics of TASEP and 
random attachment-detachment in the form of Langmuir 
kinetics (LK) displays a rich behavior including coexis¬ 
tence of low and high density regions of particles and a 
boundary condition independent phase, in the limit when 
LK competes with TASEP 0. Open TASEPs with de¬ 
fects, both point and extended, have been studied; see, 
e.g. Refs. @ which investigated the effects of the defects 
on the steady state densities and currents. In addition, 
open TASEP with a single point defect along with LK 
has been considered in Ref. 0, which finds a variety of 
phases and phase coexistences as a result of the compe¬ 
tition between the defect and LK. 

In recent studies involving asymmetric exclusion pro¬ 
cesses on closed inhomogeneous rings, the total particle 
number is held fixed by the dynamics, as expected in 
exclusion processes; see, e.g., Refs. Nonconserv¬ 

ing LK is expected to modify the steady state densities of 
pure TASEP on a closed inhomogeneous ring. TASEP on 
a perfectly homogeneous ring yields uniform steady state 
densities, due to the translational invariance of such a 
system. Evidently, introduction of the particle noncon¬ 
serving LK should still yield uniform steady state den¬ 
sities, again due to the translational invariance of the 
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system, although the actual value of the uniform steady 
state density should now depend upon LK. Nonuniform 
or inhomogeneous steady states are expected only with 
explicit breakdown of the translation invariance, e.g., by 
means of quench disorder in the hopping rates at dif¬ 
ferent sites. Studies on this model should be useful in 
various contexts, ranging from vehicular/pedestrial traf¬ 
fic to ribosome translocations along mRNA, apart from 
theoretical interests. For example, consider pedestrian or 
vehicular movement along a circular track with bottle¬ 
necks/constrictions, where overtaking is prohibited and 
pedestrians or vehicles can either leave and join the cir¬ 
cular track (say, through side roads) randomly (§], or, for 
instance, consider the motion of ribosomes along closed 
mRNA loops with defects where the ribosomes can at¬ 
tach/detach to the mRNA loop stochastically la¬ 
in this article, we introduce a disordered TASEP on a 
ring with LK, where the disorder is in the form of piece- 
wise discontinuous hopping rates across the two segments 
of the ring. The unidirectional hopping of the parti¬ 
cles across a slow segment yields reduced particle cur¬ 
rent. Evidently, this breaks the translational invariance. 
Hence, inhomogeneous steady state densities cannot be 
ruled out. In addition, we allow random attachment- 
detachment of the particles or LK at every site of the 
ring. Thus, the interplay of the quenched disorder in the 
hopping rate with the consequent absence of translation 
invariance and LK should determine the steady states 
of the model. For simplicity we assume equal rates for 
attachment and detachments. Our model is well-suited 
to analyze a key question of significance, viz, whether 
the steady state density profiles and the average particle 
numbers in the steady states can be controlled by the dis¬ 
order and (or) the LK. Recent studies of nonequilibrium 
steady states in TASEP on a ring with quenched disor¬ 
dered hopping rates without any LK revealed macroscop- 
ically inhomogeneous steady state densities in the form 
of a localized domain wall (LDW) for moderate average 
particle densities in the system; see, e.g., Refs. (jj. fl3|. 
Our work provides insight about how the steady states 
in the models in Refs. [SH3 are affected by particle non- 
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conservation and allows us to study competition between 
bulk LK and asymmetric exclusion processes in ring ge¬ 
ometry. To our knowledge, this has not been studied 
before. We generically find (i) nonuniform steady states 
and phase coexistences, (ii) phase transition between dif¬ 
ferent states of phase coexistences and (iii) the system is 
always half-filled in the steady state for the whole rele¬ 
vant parameter range, regardless of the detailed nature of 
the underlying steady state density profiles. The rest of 
the article is organized as follows. In Sec.QIl we construct 
our model. Then we calculate the steady state density 
profiles for an extended defect in Sec. 1III Al and for a 
point defect in Sec. IIIIBI In Sec. hyi we compare the re¬ 
sults for extended and point defects. Next, in Sec. El we 
discuss why the average density shows a fixed value for 
any choice of the phase parameters. Finally, in Sec. IVII 
we summarize and conclude. 


II. THE MODEL 

We consider an exclusion process on a closed ID inho¬ 
mogeneous ring with N sites, together with nonconserv¬ 
ing LK. The quenched inhomogeneity is introduced via 
space-dependent hopping rates. The parts with lower 
hopping rates can be viewed as defects in the system. 
Specifically, our model consists of two segments of gener¬ 
ally unequal number of sites. We call the parts as chan¬ 
nel I (CHI) with Ni sites (sites i = 1,2, and unit 

hopping rate, and channel II (CHII) with N 2 sites (sites 
i = M + 1, M + 2,..., N ) and hopping rate p < 1, where 
M < N (see Fig. [TJ) . We consider both the cases of an 
extended and a point defect separately. The size of an 
extended defect scales with the system size, such that 
even in the thermodynamic limit, it covers a finite frac¬ 
tion of the ring (i.e., finite N 2 /N in the limit N —> oo). 
In contrast, the size of a point defect does not scale with 
the system size, and hence, the size of a point defect rel¬ 
ative to the system size vanishes in the thermodynamic 
limit (i.e., N 2 /N —> 0 for N —> oo). Now consider an 
asymmetric exclusion process on the ring: A particle can 
hop in the anticlockwise direction to its neighboring site 
iff the site is empty. No two particles can occupy the 
same site and neither can a particle move backward even 
if there is a vacancy. Between CHI and CHII, particle 
exchanges are defined at the junctions A and B only (see 
Fig. ID , where dynamic rules are defined by the originat¬ 
ing site. In addition, the system executes LK, e.g., at any 
site, a particle can either attach to a vacant site from the 
surroundings or leave an occupied site at a rate oo. 

With nu and nm as the steady state number densities 
at i-tli sites of CHI and CHII, respectively and N to t the 
total number of particles in the system (N\ = IN, N 2 = 
(1 — l)N, l refers to the fraction of sites having hopping 
rate unity), 

[nil + n n (l - l)]N = N tot , (1) 

where, n/ = JAri/i, nu = Y^i n m ■ Now define the 



FIG. 1: (Color online) Particles hop with rates 1 and p < 1 
in segments CHI (red) and CHII, respectively, oo = fl/N (see 
text) denotes both the evaporation and condensation rates; A 
and B are the junctions where the two segments meet. 


mean number density for the total system, n as 

n= = n il + nn(l-l). (2) 

Due to LK, n is not a conserved quantity and cannot 
be used to characterize the steady states in the model, 
unlike [§j. Rather, l, p and oo parametrize the steady 
density profiles. In order to ensure that the total flux 
of the particles due to LK is comparable to the par¬ 
ticle current due to the hopping dynamics of TASEP 
(i.e., the total detachment-attachment events of the par¬ 
ticles due to LK should be comparable to the number 
of crossings of the junctions A and B by the particles 
in a given time interval), we introduce a scaling for the 
evaporation/condensation rates and define the total rate, 
= ooN and analyze the system for a given D ~ 0(1). 
This ensures competition between LK and TASEP; see, 
e.g. Ref. 1. Although there is no particle number con¬ 
servation either in Ref. f(5l] or here, it is important to em¬ 
phasize one important difference between the two that 
stems from the fact that our model is closed. As a result, 
there is no injection or extraction of particles at desig¬ 
nated ’’entry” or ’’exit” sites unlike in open TASEP or 
the model in Ref. [fj, where these rates are the tuning 
parameters. 


III. STEADY STATE DENSITIES 

We perform Mean field (MF) analysis of our model, 
supplemented by its extensive Monte Carlo Simulations 
(MCS). 
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A. MF analysis and MCS results for an extended 
defect 

Before we discuss the details of the MF analysis of our 
model, it is useful to recall the results from the model 
in Ref. where the steady state densities of a TASEP 
with open boundaries together with LK having equal 
attachment and detachment rates, 17 are investigated. 
Depending upon the entry (a) and exit (/3) rates and 
17, the steady state densities can be low or high densi¬ 
ties, or phase coexistences involving three or two phases. 
By varying the above control parameters, transitions be¬ 
tween the different steady states are observed. Motivated 
by these results and considering the junctions between 
the two segments as the effective entry and exit points 
of the segments, it is reasonable to expect similar behav¬ 
ior including phase coexistences and transitions between 
them in our model. Our detailed analysis as given be¬ 
low partly validate these expectations; we show that our 
model displays phase coexistences, but has no analogs of 
the low and high density phases. We find that even in 
the case of a point defect where there is effectively only 
one junction, these remain true. 

Our MF analysis is based on treating the model as 
a combination of two TASEPs - CHI and CHII, joined 
at the junctions A and B, respectively. Thus, junctions 
B and A are effective entry and exit ends of CHI. This 
consideration allows us to analyze the phases of the sys¬ 
tem in terms of the known phases of the open boundary 
LK-TASEP :6]. For convenience, we label the sites by 
a continuous variable x in the thermodynamic limit, de¬ 
fined by x = i/N, 0 < x < 1. In terms of the rescaled 
coordinate x , the lengths of CHI and CHII are l and 1 — l, 
respectively. For an extended defect here, l < 1. Without 
the LK dynamics, the steady state densities of a TASEP 
on an inhomogeneous ring may be obtained by means 
of the conservation of the total particle number and the 
particle current in the system |9Hlll|. In contrast, it is 
important to note that in the present model, due to the 
nonconserving LK dynamics, the particle current is con¬ 
served only locally , since the probability of attachment or 
detachment at a particular site vanishes as 1 /N Q. Sim¬ 
ilar to Ref. the steady state densities, ni{x), nn(x) 
follow 


(2 "'- 1) (^ i ‘ n ) = 0> (3) 

(2 n n - 1) (^- -Qj =0. (4) 

These yield nj, nn = 1/2 and ri/, nn = 17a; + Ci/Cn , 
where C/,C// are the integration constants. 

Apply now the current conservation locally at A and 
B. Ignoring possible boundary layers, this yields 

ni(l-ni) =pnn{l-nn) ( 5 ) 

separately, very close to x = 0 and x = l. Since p < 1 nec¬ 
essarily, Eq. © yields nj 1/2. Thus, either nj > nn or 


7ij < Tin very close to the junctions A and B. It is known 
that with 17 = 0, njj = 1/2 is a solution for moderate 
n, with 71/(a;) being in the form of a localized domain 
wall (LDW) Um- Finite 17 is expected to modify these 
solutions. Nonetheless, since nn = 1/2 is a solution of 
the steady state equation ©, nn = 1/2 remains a valid 
steady state solution for non-zero 17. Whether or not 
there are other solutions for nn is discussed later. With 
nn = 1/2, we obtain ni(x) at x = 0,1 by ©. 

Application of © yields 


71/ 


1±VT^ 

2 


( 6 ) 


at x = 0, l , which serve as boundary conditions on m(x). 
It is useful to compare CHI with an open TASEP. We 
identify effective entry (a e ) and exit (/3 e ) rates: a e = 
(3 e = (1 — \Jl — p )/2 < 1/2. Considering the fact that 
the hopping rate p of CHII is less than that in CHI (unit 
value), on physical grounds, we expect particles to accu¬ 
mulate behind junction A in CHI only. In other words, 
we expect ni{x = l) > ni(x = 0). These considerations 
allow us to set the boundary conditions ni(x = 0) = 
[1 — y/1 — p\/2 and ni(x = l) = [1 + y/1 — p|/2. Hence 
from Eq.©, we arrive at the three following solutions for 
71/(a;), namely 


nia(x) = fix + 


1-y/l-p 

2 


(7) 


Tlip(x) = ft{x — /) + 


i + yi-p 
2 


( 8 ) 


and 


nib = ^ ( 9 ) 

where ni a (x) and nip(x) are the linear density profiles 
satisfying the boundary conditions at the entrance (B) 
and exit (A) ends of CHI and nib represents the MC 
region. Notice that the solution /i/b = 1/2 cannot be ex¬ 
tended to the junctions A and B, else © will be violated. 

Given the physical expectation that ni(x) should not 
decrease with x , we identify two values of x, viz., x a 
and xp where the linear solutions meet with the third 
solution. Depending on these values of x a and xp, we will 
see that the system is found in various phases which are 
parametrized by p and 17. Using ni a (x a ) = nip(xp) = 

we get x a = and xp = l — = l — x a . Thus we 

find 

xp - x a = l- P -. (10) 

Hence there are three distinct possibilities, namely x a = 
xp, x a < xp, x a > xp, depending on which the system 
will be found in different phases. We now analyze each 
of the cases in detail. 
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(1) Consider x a < xp. Here we observe a three-phase 
coexistence. Near x = 0, we see a low density (LD) 
phase having density, ni x < 1/2, rising with a positive 
slope upto x = x a . For x a < x < xp, there is a maximal 
current (MC) phase with ni x = 1/2 and the current, 
Jix = 1/4 and while xp < x < l, we see a high density 
(HD) phase with nj x > 1/2. This is accompanied by an 
MC phase in CHII. Representative plot of comparisons 
of MFT and MC results are shown in Fig. [2] 

(2) When x a = xp, the maximal current region sepa¬ 
rating the two linear solutions vanishes and the density 
profile becomes an inclined straight line, matching con¬ 
tinuously the densities of the LD and HD phases, see 

Fig.H ' 

(3) As x a > xp , we can no more find the MC re¬ 
gion and instead find a density discontinuity. Solu¬ 
tions from the left and right meet at a point x w in 
the bulk of CHI in the form of a localized domain 
wall (LDW), where the left and right currents, i.e., 
Ja{x w ) and Jp(x w ) are equal. We can arrive at an 
expression for x w using the the local current conserva¬ 
tions. Here, J a {x w ) = ni a (x w ) (1 — ni a (x w )). Simi¬ 
larly, Jp(x w ) = nip(x w ) (1 - nip(x w )). The equality of 
Ja(x w ) and Jp(x w ) gives us the condition 

Tll a {x w ') +TH I p{x w ) — 1, (1-0 

since m a {x w ) ± nip(x w ). Using Eg. (fill). x w = 1/2. 
Thus, the LDW is always at the midpoint of CHI, unlike 
in Refs. mm- The fact that x w = 1 /2 may be un¬ 
derstood from the symmetrical structures of ni a (x) and 
nip(x). Notice that l/2 — m a (x = 0) = nip(x = Z) —1/2. 
Since, the upward slope of ni a {x) is same as the down¬ 
ward slope of nzp(x), which is H, equality of the cur¬ 
rents J a (x w ) = Jp{x w ) ensures that x w = 1/2. Thus, x w 
is independent of the values of Q and p. This is to be 
contrasted with Ref. [sj where D is known to affect the 
location of the LDW there. The height of the LDW is 

A/i e — n.jp^Xyj') nj a (.x w ') — \J 1 p HZ, (12) 


which depends upon Z, p and D. In our MCS studies 
we have generally chosen 1 = 1/2 without any loss of 
generality. Since Z = 1/2, the domain wall in CHI will 
always be located at x w =1/4, irrespective of the values 
of p and fh See Fig. 0] for a representative plot. The 
overall average density n may be found from (neglecting 
the boundary layers) 


n= ni(x)dx + / nu(x)dx. 


(13) 


Substituting the MF forms of ni(x) and nn(x), it is clear 
that n = 1/2, regardless of whether CHI is in its two- or 
three-phase coexistence state. Our MCS results agree 
with the MF results to a good extent. 

We now discuss why the system cannot be found in 
any other combination of the phases of an open TASEP. 
First of all, there is no possibility of only an MC phase 



FIG. 2: (Color online) Steady state density profiles showing 
three-phase coexistence in CHI for l = 0.5, N = 2000. The left 
and right solutions meet at n( x) = 1/2 in the bulk. Square 
points represent MCS results for p = 0.8,12 = 1.5; while the 
deep and shaded lines show the corresponding MFT equations 
for nic and nip. The circular points represent MCS studies 
for p = 0.55, 12 = 1.5. Clearly, the length of the MC region 
varies depending on the parameters, p and 12. 



FIG. 3: (Color online) Straight line density profile for the 
extended defect ( x a = xp,l = 0.5,1V = 2000) in CHI. The 
dotted black line represents the MFT denisty profile equation 
for 12 = 1.6, p = 0.36 and the triangular points shows the 
corresponding MCS result. The square points correspond to 
the MCS density profile for 12 = 1.6, p = 0.32. 


in CHI. This follows from the fact that if m = 1/2, Eq. [5] 
would be violated at the boundaries. In order to have an 
LDW in CHII, particles should pile up behind CHI which 
is physically unexpected since CHI has a higher hopping 
rate. Further, we argue that with an MC phase in CHII, 
CHI cannot be found in a pure LD phase. For CHI to be 
in such a phase, x a = 1/2. But it is also necessary that 
xp has to be greater than x a , otherwise CHI will show 
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FIG. 4: (Color online) Formation of an LDW in CHI for x a > 
xp with l = 0.5, N = 2000. Two different data sets have 
been plotted for p and fl. Points represent MCS results. The 
lines represent the following (all MFT equations): Broken 
rectangles(ni,a, II = 0.2, p = 0.2), Broken square(ni a , = 
0.2, p = 0.2), Solid(ni^, = 0.8, p = 0.2), Gapped (m a , Q = 
0.8, p = 0.2). As expected, the location of the domain wall 
remains fixed (continuous line at x = 0.25 represents the MFT 
DW). 


an LDW. But the maximum possible value for xp is 1/2. 
Now, x a = xp = 1/2 implies p = 1, for which the system 
becomes homogeneous. Since we necessarily have p < 1 
in our model, a pure LD phase for CHI with an MC phase 
in CHII is ruled out. Due to the particle-hole symmetry, 
we rule out a pure HD phase for CHI with CHII in its MC 
phase. Lastly, both CHI and CHII cannot be in their LD 
phases. This may be understood as follows. The general 
solutions for Eqs. © and © are either inclined lines of 
the form SI x + Cjjj {Cijj being constants) or flat (1/2). 
For LD phases, Ciji are to be determined by the density 
values at the entry sides of CHI and CHII, respectively. 
Let us assume nj (x = 0) = nn and nu(x = l) = nm. 
With these known values, n/( x = l) and nu(x = 1) 
respectively in CHI and CHII can be determined. Say, 
these values are n/ r and n// r , for the respective channels. 
Clearly, nj r > nn and n// r > nm. But n/ r is connected 
to nm at junction A by 

nir(l - nir) = pnm(l - nm) (14) 

and nii r and nn are connected at junction B by 

pn Ur {l - n IIr ) = nn(l - n n ). ( 15 ) 

Clearly, both m and m cannot be satisfied, simulta¬ 
neously. Hence this rules out the LD-LD phase for the 
system. Similar arguments rule out simultaneous HD 
phases in both channels. 


1. The phase diagram 


We now discuss the phase diagram spanned in the fl-p 
space. Consider the case when x a = xp. Using Eg. fTOl 


V 1 -p = 

(16) 

which gives the phase boundary : 


p = 1 — [f^]" 

(17) 


Thus when x a < xp, we have p > 1 — [DZi] 2 and the 
system shows three-phase coexistence of LD (nj < 1/2), 
MC (ni = 1/2) and HD (nj > 1/2) regions in CHI. For 
p < 1 — [LH] 2 , we get an LDW in CHI. For both three- 
phase and two-phase coexistences in CHI, CHII will be 
in its MC phase. The phase diagram is shown in Fig(5] 



FIG. 5: (Color online) The phase diagram for l = 0.5. The 

q2 

curve p = 1- — separates regions of two-phase and three- 

phase coexistences. 

The width W of the MC region (numerically equal to 
xp~x a > 0) in the three-phase coexistence, can be iden¬ 
tified as the order parameter for the phase transition be¬ 
tween a three-phase coexistence and a two-phase coex¬ 
istence. When the system is in three-phase coexistence, 
W > 0, where as it is zero in the two-phase coexistence. 
For a fixed value of p, as Q is increased, the system makes 
a transition from a two-phase coexistence state to a three- 
phase coexistence state following Eq. lfTTl) . Accordingly, 
W increases from 0 to l as H is increased for a fixed p. At 
the transition, W = 0. We now find how W approaches 
zero as p — » p c or D —> Q c in the three-phase coexistence; 
Pc tt c are located at the phase boundary given by El. 
Writings = p c +Sp, D = n c +i5H for small Sp > 0, Sfl > 0, 
we have 

W = l - ^=-^[1 - Sp/ 2 - 5Q] = V1 ~ P ° ( Sp/2 + Sn ), 

a Lq u Lq 

(18) 
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to the linear order in Sp, SQ. Evidently, W vanishes 
smoothly as Sp and S fl, indicating the second order na¬ 
ture of the phase transition. Thus, considering either p 
or n as the control parameter (for fixed f l or p, respec¬ 
tively) and drawing analogy with equilibrium second or¬ 
der phase transitions 0 we extract an ’’effective” order 
parameter exponent of value unity. This is to be con¬ 
trasted with the MF order parameter exponent of value 
1/2 in equilibrium critical phenomena. This difference 
is not surprising, considering that the present model is 
inherently out of equilibrium. 


B. Density profiles for a point defect 


Consider now the extreme limit with l —> 1, i.e., with 
a point defect. In this limit, the system has only one site 
where the hopping rate, p is less than 1 while for all other 
sites, the hopping rate is unity. Thus the MF analysis 
above by considering the system to be a combination of 
two TASEP channels joined at two ends no longer works 
because CHII (as defined for an extended defect) now 
contains just one site and has a vanishing length relative 
to the whole system for A" —> oo. Instead, the system is 
just one TASEP, say CHI, with a density ni(x), 0 < x < 1 
and two of its ends joined at one site having a hopping 
rate p < 1. 

Let the defect be present at x = 0 (which is same 
as x = 1). We assume that the particles are hopping 
anticlockwise as before. On physical grounds we expect 
piling up of particles (if at all) should occur behind the 
blockage site at x = 0. 

Now assume a macroscopically nonuniform steady 
state density profile, such that there is a pile up of par¬ 
ticles behind the defect at x = 0 and hence a jump in 
the density at a; = 0. Let pj and pjj be the densi¬ 
ties just to the left and right of x = 0, respectively: 
nj( 1 — e) = pi, ni(e) = pn,e —> 0 Using current con¬ 
servation at x = 0, we write 

Pi{ 1 - Pi) = PPi(± ~ Pii) = Pii{ 1 - Pii) ( 19 ) 

This yields solutions for pi and pn, viz., pi = -j-A and 
Pii = Density ni(x) satisfies the equation 

(2 m -1) - nj = 0, (20) 

yielding solutions 

ni(x) = 1/2, fix + C. (21) 

The constant of integration C is to be fixed by using 
either of the boundary conditions pi and pn. These evi¬ 
dently yield two values of C, say Cr and Cl, respectively, 
giving Cl = and Cr = A- _ Q. Therefore, the two 


solutions of ni(x) are 

ul(x) = fix -)-( 22 ) 

p+ 1 

n R {x) = tt(x - 1) H- 

p + 1 

along with the uniform solution rib = 1/2. Similar to an 
extended defect, we can compare CHI in case of point 
defect with an open TASEP and extract effective entry 
and exit rates: a p = /3 P = p/(p + 1) < 1/2 where a p and 
/ 3 P are entry and exit rates respectively. 

Since ul(x) and tir{x) depend linearly upon x, in gen¬ 
eral they should meet with the uniform solution, i.e., 
rib = 1/2 at two points say, xl and xr. The quantitative 
analysis follows the same logic as above for an extended 
defect. Accordingly, the values of xl and xr will deter¬ 
mine whether the system is in its three-phase coexistence 
state or a two-phase coexistence state. We find 

= (23) 


XR ~ ( 1/2 + 0 -^)! (24) 

The system will thus be in three-phase coexistence when 
xl < xr and we will have two inclined lines meeting the 
third solution in the bulk at Xl and xr, respectively. The 
extent of the MC phase (bulk solution) is given by 


xr - x L = 1 - 


1 ~P 
1 +p 


n. 


(25) 


Similarly, for xl > xr, the system will be found in its 
two-phase coexistence state and the two inclined solu¬ 
tions will meet in the bulk in the form of an LDW. The 
location of the LDW, x (j, may be calculated similarly as 
that for an extended defect, yielding x^, = The height 
of the LDW is density difference between niL(x) and 
rijR^x) at x = 1/2 which is given by 


A h p = n R (Xp) - n L {Xp) 


1 ~P 
1 +p 


- n. 


(26) 


Thus, with both an extended and a point defect, the 
location of the LDW is at the middle of CHI. Again as for 
an extended defect, this is a consequence of the symmetry 
in the forms of ul{x) and ur(x). For xl = xr, the extent 
of the MC phase vanishes and one obtains a straight line 
smoothly connecting the densities Ur and Ur. Our MFT 
results here are complemented by extensive MCS studies. 
Plots of ni{x) versus x for various p and Q in the steady 
states are shown in Figs. [6l [7] and [ 8 ] As for an extended 
defect, we find the steady state average density n = 1/2 
in the steady state by using the MF form of n 7 (x) as given 
in Eg. (1221) . again in agreement with the corresponding 
MCS results. Lastly, there are no pure LD, HD and MC 
phases in CHI for reasons very similar to the reasons for 
nonexistence of those phases in CHI with an extended 
defect, as discussed above. 
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FIG. 6: (Color online) Density plot for xl > xr. f 1 = 0.2, p = 
0.2. The square points and the triangular points show MCS 
results for N = 1000 and N = 2000, respectively. The black 
line corresponds to the position of DW according to MFT.The 
gapped and solid lines represent Ur and til MFT equations, 
respectively. It can be seen that the DW becomes sharper on 
increasing the system size. 


FIG. 8: (Color online) Density plot for three-phase coexis¬ 
tence for the point defect, fI = 0.9, p = 0.2, N = 1000. The 
gapped and solid lines represent MFT equations while the 
square points are obtained from MCS. 


by the condition xl = xr. This yields 



FIG. 7: (Color online) Straight line profile for the point de¬ 
fect (xl = xr,N = 1000). Even though qualitatively MFT 
and MCS agree, minor quantitative disagreements between 
the two for the point defect case can be seen. The dotted 
black line represents the MFT density profile equation for 
If = 0.6, p = 0.25 and the green line shows the corresponding 
MCS result. The red points correspond to the MCS density 
profile for Q = 0.6, p = 0.22. 


C. Phase diagram for a point defect 


P = 


1-0 
1 + 0 ' 


(27) 


The phase boundary is shown in Fig. [3] Following our 
analysis above for an extended defect, we consider the 
width of the MC phase, W p = xr — Xl > 0 as the or¬ 
der parameter for the phase transition between the three- 
and two-phase coexistences; W p is zero in the two-phase 
coexistence. Again as above, W p depends linearly on Sp p 
and dQ,p, where p = p C p + dp p , O = O cp — iJO p ; O cp and p cp 
satisfy Eq. J27D- Clearly, W p vanishes smoothly as Sp p 
and 5n pi indicating the second order nature of the phase 
transition. Again, the corresponding order parameter ex¬ 
ponent is unity, same as that for an extended defect. 



To construct the phase diagram in the f l — p plane, we 
identify the threshold line for crossover from three-phase 
to two-phase coexistence. The crossover line is obtained 


FIG. 9: (Color online) The phase diagram for a point defect. 
The curve p = (1 — fi)/(l + fi) separates regions of two-phase 
and three-phase coexistences. 









IV. COMPARISON BETWEEN THE DENSITY 
PROFILES WITH EXTENDED AND POINT 
DEFECTS 


We find that both with extended and point defects, 
the system can be either in three-phase or two-phase co¬ 
existence. As a result, the phase diagrams <J5]) and © 
respectively, for extended and point defects have simi¬ 
lar structures. However, the precise phase boundaries in 
the two cases differ significantly quantitatively. Similar 
differences are observed in plots of density profiles with 
both extended and point defects; for plots of density ver¬ 
sus x, see Figs. [10] and [Tl] for point defects and Fig. [12] 
for extended defect. The difference between MFT and 
MCS results is markedly visible in Fig. |TT] i.e., for suf¬ 
ficiently small fl. Fig. 9 shows that for certain range of 
values of the phase parameters p and f l, MFT and MCS 
yield contrasting results. For example, density profiles 
for a point defect ( p = 0.6, f1 = 0.2) are shown in Fig. ITT! 
These chosen values of p and fl are such that they lie 
in the region of the phase space where MFT and MCS 
phase boundaries do not overlap (see Fig. ©. These dis¬ 
agreements are clearly visible in Fig. [TTJ While the thick 
continuous lines (with LDW at x = 0.5) in Fig. [Tl] cor¬ 
respond to MFT solutions for the two-phase coexistence 
with p = 0.6, H = 0.2 (see Eqs. O , the MCS studies 
on the other hand, show three-phase coexistence for the 
same p and S2 (data points in Fig. [II]). The MCS results 
for N = 1000, 2000, 3000 show no systematic system size 
dependences and overlap with each other. Notice that the 
relative quantitative inaccuracies of MFT in a closed ring 
TASEP with a point defect without any LK, in compar¬ 
ison with a closed TASEP with extended defects (again 
without LK) have been known; see, e.g. Ref. |10| and 
Ref. El for studies on closed TASEPS with point and 
extended defects, respectively. It is, perhaps, not a sur¬ 
prise that for small some of theses mismatches are 
observed in our model as well. 

These differences may be understood in terms of the 
differences in the steady state currents across the defects 
in the two cases. In case of an extended defect, the cur¬ 
rent in the neighborhood of the junctions A and B is 
J e = p/4, where as it is J v = p/{\ + p) 2 very close to 
a point defect. Thus, J e < J p generically, since p < 1. 
Now, the existence of a three-phase coexistence requires 
the currents rising from low values ( J e or J p , as the case 
may be) to 1/4. The corresponding densities in CHI also 
vary linearly with slope Q from low or high values near 
the extended or point defects to reach 1/2 at the meet¬ 
ing points with the MC part of the three-phase coexis¬ 
tence, i.e., at x a and xp for an extended defect and xl 
and xr for a point defect. Since J e < J p generically, 
the corresponding values of ni{x) near a point defect is 
closer to 1/2 than they are near the junctions A and 
B with an extended defect: ni a (x = 0) < ul(x = 0) 
and nip(x = l) > ur{x = 1) for a given p. Now, for a 
given Q, the slope of the spatially varying solutions for 
nj(x) are same (H) for both extended and point defects. 



FIG. 10: (Color online) Plot of the density profile for a point 
defect displaying the mismatch between MFT and MCS pre¬ 
dictions (f! = 0.2, p = 0.4, N = 1000). The black line at 
x = 0.5 corresponds to the MFT LDW. Gapped and solid 
straight lines refer to MF equations for ur(x) and ul{x), re¬ 
spectively. Square points show MCS results. 



FIG. 11: (Color online) Mismatch between MFT and MCS 
predictions for the density profile with a point defect (fl = 0.2, 
p = 0.6, N = 1000 (square), 2000 (circle), 3000 (triangle)). 
While the MFT solutions, shown by the continuous black lines 
(with LDW at x = 0.5) for the given values of S2 and p dis¬ 
play a two-phase coexistence, MCS data points for all the 
three values of N display a three-phase coexistence. This dis¬ 
agreement is also visible in Fig. [0 where the MCS and MFT 
phase boundaries differ substantially for small fi. 


The densities ni a (x), nip(x) (for an extended defect) and 
nL{x),nn(x ) (for a point defect) must reach 1/2 in the 
bulk for the MC phase to exist. Hence, with an extended 
defect for a given p higher values of Q are required for 
the system to reach the threshold of existence for the 
MC part of the three-phase coexistence. This explains 
the quantitative differences in the phase boundaries in 
the phase diagrams © and ©. 

Notice that for the phase boundary corresponding to 
an extended defect (Fig. ©, the level of quantitative 
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FIG. 12: (Color online) Density profile for the system with 
an extended defect (1 1 = 0.2, p = 0.4, N = 2000) show¬ 
ing strong agreement between MCS and MFT predictions. 
The black vertical line corresponds to the MFT LDW. Thin 
gapped(ni a (o:)) and thick broken(ni^(a:)) lines refer to the 
MF equations. The square points show the MCS results. 


agreement between the MFT and MCS results for small- 
12 is much stronger than in the phase diagram (Fig. [9]) 
for a point defect. Similarly, density profiles for a point 
defect show much larger discrepancy between the MCS 
and the MFT results, in comparison with the same for 
an extended defect; see Figs. [lU] and ED for a point de¬ 
fect which clearly display the quantitative disagreement 
between MCS and MFT results. In fact, this disagree¬ 
ment is revealed vividly in Fig. 1111 where the density 
profile obtained from MCS describes a three-phase coex¬ 
istence, where as MFT predicts a two-phase coexistence. 
In contrast, the density profile for CHI with an extended 
defect as obtained from MCS match very well with the 
corresponding MFT prediction, see Fig. EH These fea¬ 
tures may be qualitatively understood as follows. In the 
small-12 limit, the effects of particle nonconservation is 
small, particle number is weakly conserved and hence the 
correlation effects due to (weak) conservation of particle 
number should be substantial, rendering MFT quantita¬ 
tively inaccurate. For an extended defect with 0 < l < 1, 
the total particle number in CHI should fluctuate within 
a range 1/2 ± (1 — l)/ 2, assuming l > 1 — l and the av¬ 
eraging occupation number to be 1/2 in CHI and CHII 
(borne out by our MCS and MFT). Thus, even in the 
small-12 limit, particle number fluctuations in CHI is ex¬ 
pected to be still substantial with an extended defect, 
weakening the correlation effects. This explains the good 
agreement between the MFT and MCS for extended de¬ 
fects. In contrast, for a point defect the total particle 
number fluctuations in CHI should vanish in the small-12 
limit, since l —> 1 or CHII has a vanishingly small size. 
Hence, the correlation effects should be large, causing 
larger discrepancies between the MCS and MFT results 
for a point defect in the small-12 limit. For larger 12, 
LK ensures lack of any correlation effect regardless of a 


point or an extended defect, leading to good agreements 
between MFT and MCS results for both of them. 


V. WHY n = 1/2 FOR ALL p AND 12? 

As we have shown above, for both point and extended 
defects the system is half-filled in the steady state, i.e., 
the global average density n = 1/2 for all 12 and p, disre¬ 
garding the possible boundary layers (which have vanish¬ 
ing thickness in the thermodynamic limit). This result 
may be easily obtained by using the MF form of ni(x) 
(and also of nu{x) in case of an extended defect) as given 
above for point and extended defects; our MCS results 
also validate this to a good accuracy. That n = 1/2 
generically in our model may be understood physically 
as follows. Notice that if we set p = 1 (homogeneous 
limit), the model is translationally invariant and the av¬ 
erage steady state density at every site is 1/2. Thus, 
there are equal number of particles and holes on average. 
A shift of n from 1/2 would indicate either more particles 
or more holes in the system in the steady state. However, 
even when inhomogeneity is introduced (p < 1), there is 
nothing that favors either particles or holes, since the in- 
homogeneity that acts as an inhibitor for the particle cur¬ 
rent, also acts as an inhibitor for the holes equally. Thus, 
it is expected to have equal number of particles and holes 
on average, i.e., n = 1/2 even with p < 1. Notice that 
this argument does not preclude any local excess of par¬ 
ticles or holes, since the particles and the holes move in 
the opposite directions, and hence the presence of a de¬ 
fect should lead to excess particles on one side and excess 
holes (equivalently deficit particles) on the other side of 
it. This holds for any 12, including 12 —> 0, and rules out 
a pure LD or a pure HD phase. This is in contrast with 
coupled TASEP and LK with open boundaries Q, where 
the boundary conditions explicitly favor more particles 
or holes, and consequently make pure LD or HD phases 
possible; in the limit of small-12, the boundary conditions 
dominate and the density profiles smoothly crossover to 
those of pure open TASEP. In contrast, in the present 
model, the density profiles do not smoothly crossover to 
those for a closed TASEP (with uniform density) for any 
small-12. Instead however, for 12 —> 0, the LDW heights 
(fT2l) and (126)) respectively for extended and point defects 
with LK smoothly reduce to the results of Ref. [Ti and 
Ref. dH, respectively with n= 1/2. 

VI. SUMMARY AND OUTLOOK 

In this article, we have constructed an asymmetric ex¬ 
clusion process on an inhomogeneous ring with LK. Our 
MFT and MCS results reveal that in the limit when 
TASEP along the ring competes with LK, the model dis¬ 
plays inhomogeneous steady state density profiles and 
phase transitions between different states of coexistence, 
parametrized by the scaled LK rate 12 and hopping rate 
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p < 1 at the defect site(s). The model always has a mean 
density n = 1/2 for all fl and p; as a result there are no 
pure LD or HD phases, in contrast to LK-TASEP with 
open boundaries. Our model may be extended in various 
ways. For instance, we may consider unequal attachment 
(to a) and detachment lud rates (to a/cod = K ^ 1). In 
this case, by using the arguments above, n = K/(l + K), 
see jgj and can be more or less than 1/2. Thus, far 
more complex phase behavior including spatially varying 
LD or HD phases should follow [g}. Details will be dis¬ 
cussed elsewhere. Additionally, one may introduce mul¬ 
tiple defect segments or point defects. However, unlike 
Refs. [IS HU, no delocalized domain walls (DDW) are 
expected even when various defect segments or point de¬ 
fects have same hopping rates. This is because DDWs 
in Refs. BI05 are essentially consequences of the strict 
particle number conservations in the models there, the 
latter being absent in the presence of LK. Our results 
very amply emphasize the relevance of the ring or closed 
geometry of the system in the presence of LK. The sim¬ 
plicity of our model limits direct applications of our re¬ 
sults to practical or experimental situations. Nonethe¬ 
less, our results in the context of traffic along a circu¬ 
lar track with constrictions or ribosome translocations 
along mRNA loops with defects, together with random 
attachments or detachments, generally show that the 
steady state densities should be generically inhomoge¬ 
neous regardless of the details of the defects. We hope 
experiments on ribosomes using ribosome profiling tech¬ 
niques [l5j] and numerical simulations of more detailed 
traffic models should qualitatively validate our results. 

A few technical comments are in order. First of all, 
notice that our MFT analysis is equivalent to consid¬ 


ering CHI as an open TASEP with selfconsistently ob¬ 
tained injection ( a e /a p ) and extraction (/3 e / /3 P ) rates, 
with a e = /3 e < 1/2 and a p = f3 p < 1/2, as obtained 
above. Given this analogy with an open TASEP, we can 
now compare our results with that of Ref. [6j, which has 
investigated an open TASEP with LK. With the con¬ 
ditions on a ejP and /3 e , P , our results should correspond 
to the the a = f3 < 1/2 line in the phase diagrams of 
Ref. [§i], where a and /? are the injection and extrac¬ 
tion rates in Ref. 6f|. Now notice that in Ref. |>] along 
the line a = /3 < 1/2 only two-phase and three-phase 
coexistences are possible, in agreement with our results 
here. Furthermore, as in Ref. [6|, two-phase coexistence 
is found for small a, /3, which in our model means small 
p, where as for large p, three-phase coexistence is found. 
Our MFT is based on analyzing the bulk density profile, 
neglecting the boundary layers. An alternative powerful 
theoretical approach to TASEP-like models with open 
boundaries have been formulated that makes use of the 
boundary layer itself, instead of the bulk [lj| . It will be 
interesting to extend these ideas to TASEP on a closed 
ring with LK. 
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